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Abstract. We simulate shock propagation in a wide range of core-collapsing presuper- 
novae: from compact WR stars exploding as SNe Ib/c through very extended envelopes 
of the narrow-line SNe Iln. We find that the same physical phenomenon of radiating 
, shocks can produce outbursts of X-ray radiation (with photon energy 3kT ~ 1 keV) 

■ lasting only a second in SNe Ib/c, as well as a very high flux of visual light, lasting for 
^\ ' months, in SNe Iln. 

in 

[ 1 Introduction 

\ Shock waves created at supernova explosions are observed mostly at the stage 

of supernova remnants due to their x-ray emission. For supernovae themselves 
shocks are not observed directly: normally the stage of shock break-out through 
the presupernova surface layers produces a short-lived transient of hard radia- 
tion. Yet those transients may have important observational consequences. More- 

■ over, in some supernovae the shock propagation in surrounding CSM may be 
\ decisive in producing their light on time-scales from days to years. 

A shock propagating down the profile of decreasing density should acceler- 
ate [1]. The acceleration ends when the leakage of hard photons from the shock 

■ front into outer space becomes efficient enough. We simulate shock propaga- 
tion in a variety of core-collapsing presupernova models using a code with hy- 
drodynamics coupled to multi-energy-group time-dependent radiation transfer. 
In our previous work on extended type II SNe [2], with radii of a few hun- 
dred Rq, we have found the peak effective temperature at the shock break-out, 
T e g ~ 1.5 x 10 5 K. For SN 1987A, with its presupernova radius of only ~ 50i? Q , 
we have got T c s ~6x 10 5 K [3]. Due to effects of scattering the maximum color 
temperature is a factor of 2-3 higher. 

In compact presupernovae, like SNe Ib/c, the shock can become relativistic 
[4] and is able to produce a burst of X-ray and even 7-ray radiation [5-7] . 

Supernovae of type Ib/c are also interesting for theory due to problems with 
their light curve and spectral modeling. Their understanding may serve as diag- 
nostics of mass loss from massive stars at the latest phases of evolution. 
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2 Compact presupernovae: SNe Ib/c 

Numerical modeling of shock breakout in SNe Ib/c was done previously using 
some simplifying approximations [8]. Our method, realized in the code Stella, 
allows us to get more reliable predictions for the outburst. Improvements in 
the theory done in the current work are: multi-energy-group time-dependent 
radiation transfer, and taking into account of (some) relativistic effects. 

A representative presupernova in our runs was a WR star built by the code 
kepler [9] (model 7A). Late light curves and spectra for this model were studied 
in [10]. 

The Model 7A has mass 3.199 M Q (including the mass of the collapsed core). 
Its radius prior to explosion is not strictly fixed because the outer mesh zones in 
kepler output actually model the strong stellar wind and are not in hydrostatic 
equilibrium. So, we fixed the radius by hand and we got four models with radii 
from 0.76 up to 2 R Q which are in hydrostatic equilibrium. Our results are 
summarized in the Table. 

Explosions with kepler [10] gave a maximum temperature of photons T ~ 
5 x 10 5 K at shock breakout. We have much finer mesh zoning at the edge of 
the star (down to 1 x 1O~ 12 M ) and better physics and we find much higher 
values of T. The left plot in Fig. 1 shows the difference between effective and 
color temperatures (labels 'e' and 'c', respectively, in the Table). One should 
note that the peak values of luminosity and temperature given in the table and 
on the left plot of Fig. 1 do not contain the light travel time correction. Its effect 
on the luminosity, i.e. smearing the peak on the time-scale of ~ R/c is shown 
on the right plot of Fig. 1. Models, labeled as 3.5N were computed by one of us 
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Fig. 1. Effective and color temperatures of emerging radiation (left). Shock breakout 
luminosity found by the hydrocode SNV (right) Dashed lines demonstrate the effect of 
averaging the light curve due the light travel time correction 



(D.K.N.) in 1992 with the equilibrium radiation diffusion hydrocode SNV. 
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Table 1. Parameters of shock breakouts 
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a prcsupernova mass and radius, respectively, in solar units. 
b kinetic energy at infinity in 10 51 ergs. 
c peak luminosity and temperature. 

d the width of the light curves at 1 stellar magnitude below L p . 

3 Shocks in CSM in Type Iln supernovae 

The narrow-line Type II supernovae (SNe Iln) are embedded in massive circum- 
stellar shells (wind) extending from tens of thousands solar radii (SN 1998S) to 
- 10 17 cm (SN 1988Z). One of the brightest SNIIn, SNII 1994W in NGC 4041, 
displays a spectrum dominated for ~ 200 days by low-ionization P-Cygni lines 
with widths of order 10 3 km/s, and at late times by narrower H Q in pure emis- 
sion. The lightcurve shows a plateau for ~ 120 days, after which the luminosity 
drops by ~ 3.5 magnitudes in V in only 12 days. The existence of such a drop, 
caused by circumstellar shells, was emphasized in paper [12]. The interpretation 
of spectra and light curve of SN 1994W leads to a coherent model in which the 
supernova interacts with a massive circumstellar shell ejected roughly 2 years 
prior to the explosion [11]. 

Main features of SN 1994W are well reproduced by the radiative shock prop- 
agating at relatively low speed for several months in the circumstellar shell. This 
shock forms its unique lightcurve with a plateau quite different from a classical 
SNell-P. (In the latter U BV light curves diverge at the end of V plateau, while 
here all colors converge.) The sudden drop of V-flux by ~ 3.5 magnitudes in 12 
days is also explained quite naturally. The results presented in Fig. 2 are com- 
puted for the supernova model with ejected mass 7M Q and huge radius of presu- 
pcrnova Rq — 2 x 10 4 i?© surrounded by a shell with the density p w i n d = 12/r in 
CGS units. The CS envelope has an outer cut-off radius R = 6.6 x 10 4 i? Q . The 
ejection of this dense extended CS envelope might be related to weak explosions 
occurring several years prior to the collapse of the core [13] (see also [14]). 
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Fig. 2. BVR light curves and T ph for the run sn94w58. 
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